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IMPROVEMENT OF MEAN ORBITAL ELEMENTS 
FOR VINTI’S SPHEROIDAL SATELLITE THEORY 


Harvey Walden 


ABSTRACT 

A method is proposed for improving the set of mean orbital elements for 
Vinti's spheroidal satellite theory in the absence of observational data, i.e., 
when the usual differential method of orbit improvement by a least-squares 
process is unavailable. An iterative procedure involving the first-order Taylor's 
series expansion of position and velocity residuals at epoch time is utilized. 
This process may be used to determine mean orbital elements directly from 
initial conditions or Keplerian osculating elements, thus eliminating the need 
for numerical factorization of two quartic polynomials arising in the inversion 
of the integrals of motion. Alternatively, the method may be applied to find the 
small corrections to the mean orbital elements following numerical factoriza- 
tion of the quartics, as carried through second order in the oblateness para- 
meter by Vinti. Numerical results are presented indicating both aforementioned 
applications of the method to actual satellite trajectories. 
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IMPROVEMENT OF MEAN ORBITAL ELEMENTS 
FOR VINTI'S SPHEROIDAL SATELLITE THEORY 


INTRODUCTION 

The problem of the determination of a set of constants of integration for use 
in a mathematical theory of orbital motion is ordinarily considered a wholly 
distinct problem in celestial mechanics from the calculation of motion for a 
given set of initial conditions or orbital elements.! In the spheroidal theory for 


satellite orbits developed by Vinti (References 1 and 2), a proced ure is jgr esented 
for deriving an accurate reference orbit of any drag-free satellite moving in the 
gravitational field of an axially symmetrical oblate^plan et As applie d to arti- 
ficial satellites of the Earth, the reference orbit takes into account exactly the 
effects of all zonal harmonic terms in the series expansion of the geopotential 
through the third, and, in addition, also accounts for the major portion of the 
fourth zonal harmonic. 


In the development of the spheroidal theory, the inversion of the integrals of 
motion involves the factoring of two quartic polynomials in order to determine 
the co-ordinates and velocities as functions of time. Vinti proposes (Reference 3) 

i 

that this factorization be carried out numerically by beginning with a zero-order 
solution corresponding to Keplerian two-body motion and solving by a method of 
successive approximations carried through second order in the oblateness param- 
eter. This method allows the constants of the motion originally suggested by 
Izsak (Reference 4), which permit exact factoring of both quartics, to be found 
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correct to second order. It is this set of orbital elements which is used in the 
equations of the final solution for the co-ordinates and velocities (Reference 5). 

In practice, an orbit improvement method is applied to reduce the magnitude 
of the observational residuals, that is, the numerical discrepancies between 
positions and velocities predicted by the analytical theory and those actually 
observed by instruments at the satellite tracking stations. For Vinti's spheroidal 
satellite theory, a method of orbit improvement through an analytic differential 
correction has been developed (Reference 6) which minimizes the observational 
residuals through the principle of least squares, yielding improved values for 
the orbital elements. Upon convergence of the iterated least-squares fitting 
process, a mean set W orbital elements, conditioned by the set of observational 
data, is produced for use in more accurate predictions. However, the applica- 
tion of this orbit improvement method presupposes the availability of observa- 
tional data required for the fitting procedure. 

Here we wish to propose an alternative method for improving the set of 
mean elements for Vinti's spheroidal satellite theory in the absence of observa- 
tional data. An iterative process involving a first-order Taylor's expansion of 
position and velocity residuals at epoch time is utilized. This procedure may 
be applied to find the small corrections to the mean orbital elements after the 
process of numerical factorization of the quartics has been carried through 
second order. It may also be used to determine mean orbital elements directly 
from initial conditions (or, equivalently, from Keplerian two-body orbital ele- 
ments, referred to as osculating elements), thus e liminat ing the need for numeri- 
cal factorization. This method, it might be noted, bears certain superficial 
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resemblances to an iterative method proposed (Reference 7) for an entirely un- 
related general perturbation theory for drag-free satellite orbits, but the differ- 
ences are substantial. The method to be described enjoys the advantage of being 
analytic in nature, rather than numerical, so that it is not necessary to rely upon 
a sufficient criterion for convergence for the iterative process involving partial 
derivatives which are never evaluated in practice, as is the case in the afore- 
mentioned method. 

FORMULATION OF THE PROBLEM 

The constants of the motion (i = 1, 2, ...., 6), which are the mean orbital 
elements for Vinti's spheroidal theory of satellite motion, include the following: 
the semi-major axis a, the eccentricity e, a parameter S corresponding to the 
square of the sine of the inclination of the orbital plane to the Equator in Kep- 
lerian two-body motion, a parameter /3 l corresponding to the negative of the 
time of passage through perigee in Keplerian motion, a parameter j3 2 corre- 
sponding to the argument of perigee in Keplerian motion, and a parameter /3 3 
corresponding to the right ascension of the ascending node in Keplerian motion. 
Let us assume that the initial co: iditions provided are in the form of rectangular 
inertial position components X 0 , Y 0 , Z 0 and velocity components X 0 , Y 0 , Z 0 , 
given for a particular time t Q . (If, in fact, the initial conditions provided are in 
the form of Keplerian osculating orbital elements, then these may be transformed 
to inertial position and velocity vector components with the usual two-body trans- 
formation equations.) By the use of the analytical equations of the orbital solu- 
tion, these six initial components determine the mean orbital elements, through 
second-order factorization of the quartics. The mean orbital elements are used 
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in turn to predict rectangular inertial position and velocity components for arbi- 
trary times t . The position and velocity vectors predicted analytically for the 
epoch time t 0 will be denoted X, Y , Z and X, Y , Z. These components will differ 
from the initial conditions because of inaccuracies in the numerical factorization 
process. Of course, if numerical factorization is not utilized, the discrepancies 
between the set of initial conditions and the predicted position and velocity based 
upon Keplerian osculating orbital elements will be considerably larger in mag- 
nitude. 

In any event, it is assumed that the differences in the co-ordinates and veloci- 
ties, as well as the required corrections to the orbital elements, are sufficiently 
small so that their squares and higher powers may be neglected. The residual 
differences in co-ordinates and velocities may be expressed by a truncated 
Taylor’s series expansion restricted to first powers as 


X 0 - X = 





Y 0 - Y = 



> ( 1 ) 


Z 


6 

L 


BZ 


Aq. 


J 
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Here the predicted values X , Y , Z , X , Y , and Z are viewed as functions of six 
independent variables, q 4 (i = 1,2, 6) which are to be improved by the 
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additive increments A (i = 1, 2 6). The time t 0 has been omitted as 

an independent variable because it remains constant throughout. The coefficients 
in the Taylor's expansion expressing the increment in the co-ordinates and 
velocities caused by a change in the orbital elements have the form of partial 
derivatives. These so-called differential coefficients are to be determined 
analytically from the mathematical equations of motion, as is shown in the fol- 
lowing. 

DIFFERENTIAL COEFFICIENTS 

There are a total of thirty- six differential coefficients to be determined, but 
half of these, viz., BX/Bq. , BY/Bq^ BZ/Bq i (i = 1, 2 6), have been evalu- 

ated previously (Reference 6) in the development of an analytic differential cor- 

• ■ • ■ • - 

rection. From the equations (Reference 5) for the velocities X, Y, and Z, one 
finds directly that 
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In the above, the subscript "i" ranges from one to six, corresponding to the 
orbital elements a, e, S, , /5 2 , and /3 3 , respectively. All partial derivatives 
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in Equations (2), (3), and (4) have been expressed explicitly previously (Refer- 
ence 6 ), with the exception of those of p and 77. Note that certain partial deriva- 
tives vanish, e.g., B a 3 /b/3 v Ba 3 /B/3 2 , Ba 3 /B/3 3 , Bp/B/3 3 , and 3 tj/9/3 3 . 

i t 

The equations for p and 17 are slightly modified from the versions originally 
presented (Reference 5) to the following forms. 

a e [(- 2 q,j) (p 2 - 2 b^ p + b 2 2 )] 
p = : sm E , ( 5 ) 

p 1 + C X 7} 4 


and 


Q a 2 (l + CjTj - C 2 t 7 2 ) % 

= ; CQS y . 

(p 2 + c 2 T7 2 ) 


( 6 ) 


Thus, for i = 1, 2, 3 (corresponding to a, e, S), 
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In Equation (7), the symbols S u and S 2i are Kronecker deltas, defined as fol- 
lows: 8 . . =1 when i = j, and 8 . . =0 whenever i / j. 

For i = 4,5 (corresponding to 0 lt /3 2 ), the simpler expression that follows 
applies. 
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The partial derivative of p with respect to /3, vanishes. 


In similar fashion, for i = 1,2,3, one finds that 
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while for i = 4, 5 , this reduces to 

i 
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The partial derivative of p with respect to fi 3 vanishes. 

All partial derivatives appearing in Equations (7) through (10) have been 
determined explicitly (Reference 6 ), so that the differential coefficients occurring 
in the first-order Taylor's expansion of the position and velocity residuals may 
now be evaluated. The simultaneous set (1) of six linear algebraic equations 
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admits a solution for the six unknowns, A q. , which are used to correct the 
orbital elements, as follows. 

t 

q'i = d, + Aq. , i = 1, 2, ..... 6. (11) 

These corrected orbital elements are used to predict analytically a new set of 
position and velocity components at epoch time and to re-evaluate the differen- 
tial coefficients in Equation (1). The process is then continued iteratively until 
the absolute values of the position and velocity residuals at epoch time, I X 0 - X| , 

! Y 0 — y! , . . . . , |Z 0 -z|, reach some sufficiently small predetermined 

values. 

NUMERICAL APPLICATIONS 

This iterative method of improving a set of mean elements for Vinti's 
spheroidal satellite theory without the use of observational data has been applied 
to actual satellite trajectories. The results of an application to determine mean 
IzSak orbital elements directly from initial conditions, i.e., without the use of 
numerical factorization to evaluate the elements to second order, is summarized 
in Table 1-A. The first column lists the initial conditions for a particular epoch 
of an artificial satellite orbit in the form of inertial rectangular co-ordinates 
and their time derivatives. The second column provides the osculating Keplerian 
orbital elements corresponding to the initial conditions. These classical Kep- 
lerian elements were obtained through use of the two-body equations of motion, 
considering only the central term in the potential field. The final column displays 
the converged values for the mean set of Izsak orbital elements obtained after 
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three iterations. A measure of the degree of improvement in the orbital elements 
is provided by Table 1-B, which displays the residuals in position and velocity 
components, using the osculating Keplerian elements initially and then the mean 
iterated Izsak elements obtained upon convergence. 

The results of an application to improve the mean Izsak orbital elements fol- 
lowing the use of numerical factorization of the quartics through second order is 
summarized in Table 2- A. The initial conditions for an artificial satellite orbit 
are identical to those presented previously. The second column now provides the 
mean Izsak orbital elements determined by numerical factorization carried through 
second order in the Earth's oblateness parameter, and the final column gives the 
converged values for the mean Izsak orbital elements obtained after two iterations. 
Table 2-B shows the residuals in position and velocity components with use of the 
two respective sets of mean Izsak orbital elements. 

It is seen that the final iterated set of mean Izsak orbital elements is virtually 
independent of whether the process of numerical factorization of the quartics is 
utilized. Hence the iterative method of determining mean elements for Vinti's 
spheroidal satellite theory presented in this paper may be used as a valid alter- 
native to the procedure of numerical factorization. However, if the latter proce- 
dure is applied, then substantial corrections to the orbital elements are not to be 
expected through use of the iterative method, inasmuch as any such corrections 
account for errors of third or higher orders. Without use of numerical factoriza- 
tion, the decrease in size of the position and velocity residuals through application 
of the iterative method is impressive, but even when the factorization process is 
employed, most residuals experience noteworthy decreases after only two itera- 
tions of this method. 
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TABLE 1-A 


Determination of Mean Izsak Orbital Elements Directly from Initial Conditions 


Si/ 

Initial Conditions — 

Osculating . 

Keplerian Elements- 

Mean Iterated , 
Izsak Elements — 

x = + 0.86773200 

a = 1.7478170 

a = 1.7461661 

y = + 1.0052368 

e = 0.23597612 

e = 0.23553637 

z = - 0.14256217 

sin 2 i = 0.52617614 

S = 0.52593981 

x = - 0.54766917 

-r = 0.056312786 

P x = 0.052921090 

y = + 0.38465985 

o>= 3.2487135 

yS 2 = 3.2503452 

z = - 0.69095995 

fi = -2.3849592 

P 3 =-2.3844409 


q/ 

— The position components x, y, z are in Earth equatorial radii (e.r.), and 
the velocity components x, y , z are in e.r./c.u.t., where 1 canonical unit of time 
(c.u.t.) is equal to 806.823 seconds. The semi-major axis a is in e.r., the time 
of perigee passage r and p x are in c.u.t., and the argument of perigee <^>, the 
right ascension of the ascending node 0 , /3 2 , and fi 3 are in radians. 


TABLE 1-B 

Magnitude of Residuals in Position and Velocity Components 

Residual — ^ 

Initially 

(Keplerian Elements) 

Upon Convergence 
(after 3 iterations) 

|Ax| 

1722.44 

0.29 

|Ay| 

2133.78 

0.38 

| A z | 

5425.18 

0.05 

1 Ax | 

95.2453 

0.0118 

|Ay| 

59.9060 

0.0471 

1 A z | 

85.3974 

0.0177 


g/ 

— The position residuals are in meters, and the velocity residuals in centi- 
meters per second. 
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TABLE 2- A 


Improvement of Mean Izsak Orbital Elements After Use of Second-Order 
Numerical Factorization 

, Numerically Factored , Mean Iterated , 
Initial Conditions — Mean Izsak Elements — Izsak Elements — 


+ 0.86773200 

a = 1.7461655 

a = 

1.7461661 

+ 1.0052368 

e = 0.23553626 

e = 

0.23553640 

- 0.14256217 

S = 0 .52594027 

S = 

0.52593981 

- 0.54766917 

/3 X = 0.052921097 

fii = 

0.052921090 

+ 0.38465985 

P 2 = 3.2503450 

02 = 

3.2503453 

- 0.69095995 

yS 3 =-2.3844408 

^3 = 

-2.3844408 


- The position components x, y, z are in Earth equatorial radii (e.r.), and 
the velocity components x, y, z are in e.r ./c.u.t., where 1 canonical unit of time 
(c.u.t.) is equal to 806.823 seconds. The semi-major axis a is in e.r., the param- 
eter /3, is in c.u.t., and the parameters and /3 3 are in radians. 

TABLE 2-B 

Magnitude of Residuals in Position and Velocity Components 


/ Initially Upon Convergence 

Residual — (Factored Elements) (after 2 iterations) 


| A X | 

0.0475 

0.3802 

|Ay| 

1.9008 

0.2851 

| A z ' 

1.4375 

0.0233 

| Ax | 

0.2238 

0.0118 

|Ay| 

0.0294 

0.0353 

|A z | 

0.2827 

0.0353 

a/ 

The position residuals are in meters, and the velocity residuals in centi- 


meters per second. 
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